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LIFETIME  PREDICTION  FOR  SATELLITES  IN  LOW-INCLINATION  TRANSFER  ORBITS 

by 

D.  G.  King-Hele 

SUMMARY 

a  geostationary  satellite  usually  reaches  its  final  circular  orbit 
via  a  transfer  orbit  having  apogee  height  near  36000  km  and  perigee  height 
less  than  600  km.  The  population  of  discarded  rockets  left  in  these 
transfer  orbits  is  steadily  increasing  (there  are  now  more  than  a  hundred) , 
and  their  likely  lifetimes  are  usually  assessed  after  lengthy  numerical 
integrations  to  evaluate  lunisolar  perturbations.  This  paper  gives  a 
simple  analytical  method  for  predicting  the  lifetimes  when  the  orbital 
inclination  is  between  20°  and  30°,  as  often  happens.  Several  approxima¬ 
tions  are  made,  but  the  errors  should  not  exceed  those  from  other  sources. 

Departmental  Reference:  Space  607 

Copyright 

Controller  HMSO  London 
1981 


2 


LIST  OF  CONTENTS 


1  INTRODUCTION 

2  LIFETIMES  IN  THE  ABSENCE  OF  LUNISOLAR  PERTURBATIONS 

3  CHANGES  IN  PERIGEE  HEIGHT  DUE  TO  LUNISOLAR  PERTURBATIONS 

4  WEIGHTED  MEAN  PERIGEE  HEIGHT  FOR  USE  IN  LIFETIME  PREDICTION 

5  EFFECT  OF  CHANGING  THE  INCLINATION 

6  TESTING  THE  THEORY  AGAINST  NUMERICAL  INTEGRATION 

7  ACCURACY 

8  AN  ALTERNATIVE  PROCEDURE 

9  CONCLUSION 
References 
Illustrations 

Report  documentation  page 


Page 

3 

3 

4 
6 
7 

9 
9 

1C 

10 
1  1 


Figures  1-6 
inside  back  cover 


V  Iim -if -mi-  ‘  "  '  .  .  . 


81  1  1 


3 


1  INTRODUCTION 


Many  geostationary  satellites  proceed  to  their  final  circular  equatorial  orbits  at 
a  height  near  36000  km  by  way  of  transfer  orbits  having  perigee  heights  near  200  km, 
apogee  heights  near  36000  km  and  inclinations  of  20°-30°.  A  discarded  rocket  usually 
remains  in  this  transfer  orbit,  and  its  lifetime  is  often  difficult  to  estimate  because 
lunisolar  perturbations  substantially  alter  the  perigee  height  as  time  passes.  If  the 
perigee  height  is  pushed  downwards  initially,  the  rocket  may  decay  within  a  year;  if  the 
perigee  height  is  pushed  upwards  initially,  the  lifetime  may  be  ten  years  or  more. 

To  estimate  the  likely  lifetimes  of  satellites  in  such  orbits,  it  is  usual  to  under¬ 
take  extensive  numerical  integrations,  which  not  only  eat  up  computer  time  but  may  also 
suffer  inaccuracy  because  the  initial  orbits  are  often  in  error  by  more  than  10  km  in 
perigee  height.  In  this  paper  approximate  analytical  methods  are  devised  for  estimating 
lifetime,  bearing  in  mind  the  likely  10  km  error-level  of  the  initial  orbit. 

Most  of  the  objects  to  which  the  method  applies  are  upper-stage  rockets  of 
satellites  in  series  such  as  Intelsat,  Fleetsatcom,  Westar,  Kiku,  GOES  and  SBS.  About 
forty  such  objects  are  currently  in  orbit*,  and  eight  were  launched  during  the  first 
nine  months  of  1981. 

2  LIFETIMES  IN  THE  ABSENCE  OF  LUNISOLAR  PERTURBATIONS 

A  typical  transfer  orbit  has  a  perigee  height  of  about  200  km  and  apogee  height  near 
35800  km.  This  means  that  the  semi  major  axis  a  is  near  24400  km,  the  orbital  period  T 
is  near  0.438  day  and  the  eccentricity  e  is  near  0.73. 


In  the  absence  of  lunisolar  perturbations,  the  lifetime  L  of  a  satellite  moving  in 
an  orbit  of  high  eccentricity  subject  to  air  drag  is  given  by  equation  (4.132)  of  Ref  2  as 


(1  -  e) Tf ( e)  (\  +  e\* 
4p^S  \2itHae  j 


(0 


where  f(e)  is  a  function  of  e  given  in  Table  4.2  of  Ref  2  as  1.50  for  e=0.73;  and 
H  are  the  density  and  density  scale  height  at  perigee;  and  S  is  a  drag  parameter  given 
by  6 = FSCp/m  .  Here  F  is  a  factor  allowing  for  atmospheric  rotation,  which  has  a  value 
near  0.92  for  low-inclination  transfer  orbits,  Cp  is  the  drag  coefficient,  which  may  be 
taken  as  2.2  for  heights  near  200  km,  and  m/S  is  the  mass/area  ratio  of  the  satellite. 


Inserting  the  numerical  values  for  the  transfer  orbit  into  equation  (1),  we  obtain 


L  = 


23.6 

109p 

P 


^ ioosy 


years. 


(2) 


with  m  in  kg,  S  in  m 


in  kg/m  and  H  in  km.  A  typical  value  of  m/S  is 


100  kg/tn  and  we  may  take  this  as  a  standard  value  to  give  a  standard  lifetime  Lr 


L 


0 


23.6 

I09p 

P 


years. 


(3) 
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With  values  of  p  and  H  from  the  COSPAR  International  Reference  Atmosphere 
3  P  ... 

1972  for  exospheric  temperatures  of  700  K  (low  solar  activity)  and  1000  K  (high 

solar  activity),  the  variation  of  Lq  with  perigee  height  is  as  given  in  Fig  J.  High 

solar  activity  is  not  usually  maintained  for  more  than  about  three  years,  and  low  solar 

activity  does  not  last  for  more  than  about  six  years  in  each  sunspot  cycle.  So  if  the 

calculated  lifetime  exceeds  six  years,  the  intermediate  curve  labelled  'medium  solar 

activity'  should  be  used. 

Although  100  kg/m2  is  a  useful  standard  value  for  m/S  ,  the  real  value  may  be 

considerably  different;  the  lifetime  Lq  then  needs  to  be  multiplied  by  m/lOOS  .  In 

calculating  m/S  for  a  particular  rocket,  m  should  of  course  be  taken  as  the  mass 

after  all  the  propellant  has  been  consumed,  and  the  mean  cross-sectional  area  S  is 

usually  slightly  less  than  the  side  area  id  of  the  rocket,  if  it  is  a  cylinder  of 

2 

length  l  and  diameter  d  .  If  l/d  >  2  ,  the  mean  value  of  S  may  be  taken  as 


£d(0.82  +  d/4  i)  . 

3  CHANGES  IN  PERIGEE  HEIGHT  DUE  TO  LUNISOLAR  PERTURBATIONS 

In  the  absence  of  air  drag,  the  change  in  eccentricity  due  to  lunisolar  perturba¬ 
tions  is  given  explicitly  by  G.E.  Cook  as  the  sum  of  !5  terms  in  equation  (55)  of  Ref  4. 
If  we  insert  numerical  values  for  the  orbital  inclination  of  the  satellite  (25°)  and  the 
obliquity  of  the  ecliptic  (23.4°),  and  also  take  the  lunar  inclination  as  23.4°,  equation 
(55)  may  be  written,  on  omitting  terms  less  than  21  of  the  main  term,  as 


15K 

4n 


e(  1  -  e2)  ^|o.835  sin  2(0  -  -  u^  +  m)  -  0. 147  sin(0  -  0^  +  2m)  +  0.072  sin  2(12-12^  +  0)) 

-  0.153  sin(2uj  -  12  +  12^  -  2m)  +  0.068  sin  2m j  ,  (4) 

where  K  =  0.971  deg2/day2  for  the  Sun,  K  =  2.132  deg2/day2  for  the  Moon  and  n  =  1 /T  = 

2.28  rev/day  =  821  deg/day.  In  equation  (4),  0  and  12^  are  the  right  ascensions  of  the 
ascending  node  of  the  satellite  and  of  the  disturbing  body  respectively,  m  is  the 
satellite's  argument  of  perigee,  and  u^  is  the  longitude  of  the  disturbing  body. 

The  changes  in  eccentricity  due  to  lunisolar  perturbations  are  always  small,  so  that 

e  can  be  taken  constant  on  the  right-hand  side  of  equation  (4).  This  assumption  holds 

even  in  the  presence  of  drag,  because  the  lifetime  will  already  be  decided  by  the  time 

2  i 

when  e  decreases  to  0.67,  and  since  0.497  <e(l-e  )2  <0.500  for  0.67  <  e  <  0.73  , 

2  1 

we  may  conveniently  take  e(l-e  )2  =  0.5  . 

Consider  first  the  solar  perturbation.  With  the  numerical  values  already  quoted  for 
K,  n,  etc,  the  term  outside  the  curly  brackets  in  equation  (4)  becomes  0.00222  deg/day  ■ 
0.0141  rad/yr.  Also  u^  is  replaced  by  the  solar  longitude  1  ,  and  12^  =  0  .  Thus 

e  .  =  0.014l<0.835  sin  2(12  +  m-X)  -  0.147  sin(12+2m)  +  0.072  sin  2(12  +  m) 

solar  l 


0.153  sin(2X  -  12  -  2m)  +  0.068  sin  2mj 


(5) 


■■■i  in  rifjifiT*- 
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For  a  transfer  orbit  at  an  inclination  near  25°,  0  =*  -0.4  deg/day  -2.5  rad/yr,  and 
u  =*  0.7  deg/day  =*  4.4  rad/yr.  Also  X  =  6.?  rad/yr.  So,  preparing  for  integration  by 
inserting  the  rates  of  change  divided  by  their  numerical  values,  we  have 

4  ,  «  -  0. 001 4  1  ]0. 949  |2(0  +  <L  -  X)}  sin2(n  +  w-X)  +  0. 233(0 + 2w)  sin(0  +  2u) 

solar  L  1  ' 

-  0.  189  )2(0  ♦<!>)}  sin  2(0  +  w)  +  0.243(2X  -  0  -  2uiy  sin(2X  -  ft  -  2u) 

-  0. 077( 2(L)  sin  2ul 


Integrating,  we  find  the  variation  of  e  with  time  t  due  to  solar  perturbations  as 


(e_Ko) 

sole 


=  0.00141  0.949  cos  2(0  +  uj-X)  +  0.233  cos(fi+2w)  -  0.189  cos  2(0  +  u) 


+  0.243  cos  (SI  +  2ui  -  2X)  -  0.077  cos  2oj 


where  suffix  0  denotes  initial  values. 

Now  consider  the  lunar  perturbations.  On  inserting  numerical  values,  the  term  out¬ 
side  the  square  brackets  in  (4)  becomes  0.0310  rad/yr.  Since  u^  =  84  rad/yr  for  the  Moon, 
the  first  term  in  equation  (4)  has  a  frequency  of  2  weeks  and  its  amplitude  is  less  than 
10%  of  the  main  solar  variation  in  (7).  So  it  may  be  ignored  in  our  long-term  analysis. 
The  fourth  term  in  (4)  can  u1  so  be  ignored,  being  of  short  period  and  much  smaller 
amplitude.  On  integrating  the  three  remaining  terms  in  (4),  we  find,  in  analogy  with  (7), 


(e-eQ)  =  0.0031ojo.233  cos(0  -  0rf  +  2w) 

'  7lunar  L 


-  0.189  cos  2(0-0,  +  u>)  -  0.077  cos  2m  .  (8) 

d  J0 

As  ft,  has  a  mean  value  of  2ero  over  18.4  years,  and  Ift,  I  never  exceeds  13°,  we  may 
a  a 

conveniently  take  ft ,  = 0  in  (8).  Adding  (7)  and  (8)  gives  the  total  lunisolar  perturba- 
a 


-  [■ 


|o. 00134  cos  2(0  +  oj  -  X)  +  0.00105  cos(ft~2a ,) 

-  0.00085  cos  2(0  +  u>)  +  0.00034  cos(0  +  2w  -  2X)  -  0.00035  cos  2uil  .  (9) 

J0 


On  multiplying  by  a  (=  24400  km)  to  convert  from  eccentricity  to  perigee  height  , 
and  rounding  off  to  the  nearest  km,  we  have: 


yP  "  yPo 


133  cos  2(ft  +  oj  -  A)  +  26  cos\ft+2to) 


-  21  cos  2(ft  +  u>)  +  8  cos  (ft  +  2cjo  -  2  A)  -  9  cos  2to  I  km. 
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From  equation  (10),  the  mean  perigee  height  yp  over  several  cycles  is  given  by 

yp  =  ypQ  +  33  cos  2(0q  +  u>0-Aq)  +  26  cos(0q  +  2wq)  -  21  cos  2(0^  +  ^) 

+  8  cos(0q  +  2o)q  -  2Aq)  -  9  cos  2o)q  km.  (II) 


As  a  first  approximation  we  may  ignore  the  last  two  terms,  of  amplitude  8  and  9  km,  and 
also  the  effects  of  odd  zonal  harmonics"*,  which  produce  an  oscillation  with  an  amplitude 
of  about  3  km.  Then 

yp  =  yp0  +  %1  +  %2  +  %3  say*  (,2) 

where  Ay  =  33  cos  2(S2.  +  un  -  A.)  ,  Ay  .  =  26  cos(fi_ + 2w.)  and  Ay  .  =  -21  cos  2(0.  +  to.) 

pi  u  u  u  pz  U  U  pj  u  u 

are  given  in  Fig  2.  Thus,  if  the  initial  orbital  elements  of  the  satellite  are  known,  we 

can  calculate  y  .  =  a„(l  -  e.)  -  6378  km,  and  then  calculate  y  from  (12)  or,  more 

pO  0  0  p 

accurately,  from  (11). 

4  WEIGHTED  MEAN  PERIGEE  HEIGHT  FOR  USE  IN  LIFETIME  PREDICTION 

The  arithmetic  mean  perigee  height  y  is  not  the  correct  mean  value  to  use  in 

P 

calculating  lifetime,  because  the  variation  of  density  with  height  is  strongly  exponential, 
and  not  linear.  For  estimating  lifetime  we  need  to  use  the  height  y*  at  which  the 

perigee  density  pp  has  its  mean  value  p*  ,  such  that 


P 


0 


(13) 


over  a  complete  period  P  of  the  oscillation  in  perigee  height. 

approximate  the  variation  (10)  as  a  simple  sinusoid  of  amplitude 

and  frequency  f  ,  say,  we  have  y  =  y  +  47(1  -cos  ft)  ,  where 

p  m 

perigee  height.  Since  the  perigee  density  is  greatest  at  y  =  y 

pm 

height  should  be  expressed  in  a  form  which  is  correct  at  and  near 


If  we  now  further 

/332  +  262  +  21 2  =  47  km 

y  is  the  minimum 
m 

,  its  variation  with 

y  =  y  ,  that  is 
p  m 


P 


P 


P 


m 


1  -  cos 


(14) 


where  H^  is  the  density  scale  height  at  height  ym  .  On  writing  ft  =  t  and  integrat¬ 
ing  from  0  to  2it,  equations  (13)  and  (14)  give: 
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r 


where  I^(z)  is  the  Vessel  function  of  the  first  kind  and  imaginary  argument  of  degree 
tero  and  argument  z  .  Solving  (16)  for  y*  gives 


y*  =  y  +  47  -  H  in^I  (47/H  )1 
m  m  (  0  m  J 


(17) 


Since  y  +47 
'  m 


ue  may  rewrite  (17)  as 


i 


i 

i 


y* 


(18) 


where  Ay  .  =  H  inil  (47/H  ) 
p4  m  |  0  m 

Thus  the  procedure  for  lifetime  prediction  is  to  evaluate  y^  from  equation  (11)  by 

calculation,  or  (less  accurately)  from  (12)  and  Fig  2,  and  then  to  find  y*  from  (18)  by 

reading  off  Ay  ,  from  Fig  3.  The  lifetime  of  the  satellite  can  then  be  obtained  from 
P+ 

Fig  I,  taking  the  perigee  height  in  Fig  1  to  be  y*  . 

A  further  simplification  of  equation  (I!)  is  possible  if  is  near  0  or  180°,  as 

often  happens  because  the  apogee  of  the  transfer  orbit  is  usually  required  to  be  near  the 
equator.  If  u)q  =  0  or  1 80° ,  equation  (II)  reduces  to 


and  is  plotted  against  y  in  Fig  3. 


>’p  =  y  ^  +  33  cos  2(i’0  -  Xg)  +  26  cos  0^  -  21  cos  2Q ^  +  8  cos(flp-2Xp)  -  9  km.  (19) 

Although  the  approximate  method  is  intended  primarily  for  satellites  with  perigee 
height  near  200  km,  there  are  some  objects  (such  as  Comstar  rockets)  in  transfer  orbits 
with  a  much  higher  perigee,  near  600  km.  The  method  can  also  be  applied  to  these  objects, 
since  the  changes  in  the  numerical  factors  are  not  significant.  So,  to  complete  the 
picture.  Fig  4  gives  the  variation  of  lifetime  with  perigee  height  for  transfer  orbits 
with  perigee  heights  between  300  and  700  km.  The  curved  labelled  T^ =  900  K  would  be 
applicable  if  the  average  future  solar  activity  remains  at  the  levels  prevalent  during 
the  20th  century,  but  the  prediction  of  solar  activity  for  up  to  100000  years  ahead  is  a 
matter  of  some  uncertainty. 

5  EFFECT  OF  CHANGING  THE  INCLINATION 

The  analysis  has  been  pursued  on  the  assumption  that  the  orbital  inclination  i  is 
23°,  but  there  will  be  a  range  of  inclination  over  which  the  results  are  likely  to  be 
within  the  required  limits  of  accuracy.  The  validity  of  the  method  at  other  inclinations 
can  be  assessed  by  looking  at  the  three  main  components  of  the  variation  in  perigee 
height . 


The  variation  of  perigee  height 
equat ion  (10)  is 


given  by  subtracting  equation  (1!)  from 


y  -  y 

p  p 


-  33  cos  2(ft  +  ui-A)  -  26  cos(fi  +  2w) 

+  21  cos  2 ( S i  +  o))  -  8  cos(fi  +  2o)  -  23)  +  9  cos  2io  km.  (20) 
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The  first  term  on  the  right-hand  side,  of  amplitude  33  km,  has  a  period  of  about  9  months; 
the  second  term,  of  amplitude  26  km,  has  a  period  of  about  12  months;  and  the  third  term, 
of  amplitude  21  km,  has  a  period  of  about  20  months.  The  two  smaller  terms  have  periods 
of  12  and  9  months  respectively.  If  air  drag  seriously  affects  the  orbit,  these  regular 
sinusoidal  oscillations  will  of  course  gradually  be  modified,  but  if  the  lifetime  is 

5  years  or  more,  the  first  few  cycles  will  be  little  affected  by  drag  and  equation  (20) 

.  .  o 

should  be  realistic,  for  l  =  25 

A  change  in  inclination  from  25°  affects  both  the  periods  and  amplitudes  of  the 
oscillations.  Increasing  the  inclination  reduces  the  period  and  amplitude  of  the  first 
term,  and  increases  those  of  the  second  and  third  terms.  If  the  inclination  is  increased 
to  30°,  the  three  amplitudes  (33,  2b  and  21  km)  become  29,  35  and  25  km  respectively, 
while  the  periods  change  from  about  4,  12  and  20  months  to  about  8,  14  and  25  months. 

If  the  inclination  is  reduced  to  20°,  the  amplitudes  become  37,  19  and  18  km,  and  the 
periods  become  9,  11  and  17  months.  Thus,  none  of  the  amplitudes  changes  by  more  than 
9  km,  and  the  root  sum  of  the  squares  of  the  amplitudes  does  not  change  by  more  than  5  km. 
So  it  is  apparent  that  equation  (20)  can  be  used  for  any  inclination  i  between  20°  and 
30°:  the  errors  will  be  somewhat  greater  when  i  departs  from  25°,  but  should  still  not 
exceed  10  km  (sd),  and  can  be  reduced  by  using  revised  amplitudes  in  accord  with  those 
quoted  above.  Also,  of  course,  and  u  change  at  different  rates  if  i  departs 

from  25°. 

If  the  inclination  is  not  too  far  outside  the  range  20°-30°,  it  is  still  possible  to 
use  an  equation  of  the  form  (20)  but  with  the  amplitudes  adjusted  to  suit  the  inclination. 
If  we  rewrite  the  three  main  terms  in  equation  (20)  as: 

y^-y^  =  -  A j  cos  2 (h  +  w  -  A)  -  A  cos (0  +  2u>)  +  A^  cos  2(B+u)  ,  (21) 

the  amplitudes  A( ,  A^  and  A^  vary  with  inclination  in  the  manner  shown  in  Fig  5.  The 
lower  limit  of  inclination  in  Fig  5  is  extended  to  10°  so  as  to  cover  the  Ariane  third 
stage  1 98 1  —  5 7 B ,  and  30°  is  a  convenient  upper  limit  because  there  are  as  yet  no  rockets 
in  transfer  orbits  at  inclinations  between  30°  and  40°. 

Fig  5  can  of  course  also  be  used  in  the  procedure  for  lifetime  prediction.  If  the 
inclination  is  outside  the  range  20°-30°,  the  mean  perigee  height  should  be  calculated 
using  equation  (11)  with  values  of  Aj ,  A^  and  A^  from  Fig  5  to  replace  the  values  33, 
26  and  21  km  given  in  equation  (II).  (The  same  procedure  could  also  be  used  for  inclina¬ 
tions  between  20°  and  30°;  but,  as  already  discussed  earlier  in  this  section,  the  improve¬ 
ment  in  accuracy  would  only  be  marginal.)  Fig  5  also  shows  that  </a^  +  A^  +  A^  remains 
between  45  and  52  km  for  inclinations  between  10°  and  30°,  so  there  is  no  need  to  modify 
the  analysis  in  section  4.  When  i  changes  to  10°  the  amplitudes  of  the  fourth  and 
fifth  terms  in  equation  (20)  change  from  8  and  9  km  to  5  and  1  km  respectively,  so  they 
can  be  ignored. 

The  method  developed  here  fails  when  the  orbit  is  near-resonant,  as  happens  at 

,  i  o  *  •  ,  o  •  • 

l  =*  46  when  ft  +  u>  —  0  and  A^  becomes  very  large;  or  at  l  —  56  when  ft  +  2w  -  0 

,  0  • 

and  becomes  very  large.  Another  resonance  occurs  at  1  =  63.4  when  to  becomes 
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very  small  and  the  amplitude  of  the  term  in  cos  2w  in  equation  (20)  therefore  becomes 
very  large:  this  is  the  situation  for  the  Molniya  satellites*’.  Because  of  these  three 
resonances,  the  approach  adopted  here  only  succeeds  for  inclinations  up  to  about  35°. 

6  TESTING  THE  THEORY  AGAINST  NUMERICAL  INTEGRATION 

To  conduct  a  test,  we  need  a  satellite  with  a  lifetime  of  at  least  5  years  and 
inclination  near  23°.  The  example  chosen  is  197 8—7 1C,  the  third  stage  of  ESA-Geos  2, 
which  was  launched  in  July  1978  and  was  left  in  a  transfer  orbit  of  inclination  25. 9°, 
eccentricity  0.729  and  semi  major  axis  24300  km:  these  values  are  close  to  those  assumed 
in  developing  the  approximate  method. 

Fig  6  shows,  as  an  unbroken  line,  the  variation  of  in  1979-80  as  given  by 

numerical  integration  with  the  program  PROD'7,  with  5-day  integration  intervals.  The  small 
lunar  perturbations,  of  amplitude  about  4  km  and  period  2  weeks,  have  been  smoothed  out. 
(For  pictorial  examples  of  these  oscillations,  see  Fig  4  of  Ref  8,  or  Ref  9.)  The  broken 
line  in  Fig  6  shows  the  variation  of  y^  given  by 

y  =  240  -  33  cos  2 (0  +  o-i)  -  26  cos(0  +  2ui)  +  21  cos  2(0  +  u)  km,  (22) 

P 

with  the  values  of  fi  and  w  taken  from  PROD  (though  a  linear  variation  of  0  and  ui 
would  have  been  adequate). 

The  forms  of  the  two  curves  in  Fig  6  are  very  similar.  The  difference  between  them 
has  a  mean  numerical  value  of  8  km  and  a  maximum  of  18  km.  This  is  as  would  be  expected 
when  terms  of  amplitude  8  and  9  km  have  been  ignored. 

7  ACCURACY 

There  are  two  main  orbital  errors  in  the  use  of  the  approximate  method.  First, 
there  is  the  error  due  to  the  approximation,  which  will  be  largest  when  the  initial  point 
happens  to  be  on  a  relatively  ill-fitting  region,  such  as  August  1979  or  February  1980  in 
Fig  6.  The  error  from  this  source  is  probably  about  10  km  (sd),  and  can  be  reduced  if 
necessary  by  adding  the  two  extra  terms  on  the  right-hard  side  of  equation  (11).  Second, 
and  probably  worse,  there  are  the  errors  of  the  initial  orbit.  The  source  of  the  initial 
orbit  is  sometimes  a  transfer  orbit  of  the  joined  rocket-and-satel 1 ite :  this  is  accurately 
known  to  the  launch  authority,  but  it  may  be  unreliable  if  the  kick  at  separation  diverts 
the  discarded  rocket  into  a  slightly  different  orbit.  Alternatively,  the  orbit  may  be 
derived  from  subsequent  radar  observations;  but  unless  the  satellite  is  observed  near 
perigee,  it  is  extremely  difficult  with  orbits  of  e  =■  0.7  to  determine  perigee  height 
accurately  from  normal  radar  observations.  Bias  errors  of  20  km  or  more  can  easily  occur. 

The  errors  due  to  both  these  sources  can  be  checked  and  reduced  by  recalculating  the 
lifetime  with  an  'initial'  orbit  a  few  months  later,  unless  the  object  has  been  lost  - 
which  is  often  a  signal  of  a  rapid  decay  rate. 

Other  errors  arise  in  estimating  the  mass/area  ratio  of  the  satellite  and  in  pre¬ 
dicting  the  average  solar  activity  during  the  lifetime.  These  would  both  be  expected  to 
give  errors  of  10-20%  in  lifetime,  unless  the  size,  shape  and  weight  of  the  rocket  are 


...i. 


unknown,  when  the  errors  would  be  larger.  It  would  then  he  appropriate  to  use  the 
observed  rate  of  decay  during  the  early  days  in  orbit  to  estimate  the  mass/area  ratio. 

8  AN  ALTERNATIVE  PROCEDl'RK 

If  a  satellite  in  transfer  orbit  is  successfully  tracked  for  several  cycles  of  the 
perigee  height  oscillation,  for  about  5  years,  the  most  accurate  method  of  calculating 
lifetime  is  to  use  the  mean  observed  orbital  decay  rate  over  the  5  years  in  conjunction 
with  the  normal  graphical-analytical  methods  for  predicting  lifetime*^.  This  procedure 
avoids  the  effects  of  errors  in  the  initial  orbit  and  mass/area  ratio,  though  the  errors 
due  to  uncertainty  in  future  solar  activity  are  increased.  The  procedure  also  suffers 
from  two  limitations  which  do  not  affect  the  method  developed  in  this  paper:  the  satellite 
may  be  lost,  or  may  decay,  before  the  procedure  can  legitimately  be  applied. 

4  CONCLUSION 

An  approximate  analytical  method  has  been  developed  for  predicting  the  lifetimes  of 
objects  in  transfer  orbits  with  perigee  height  between  120  and  300  km,  apogee  height  near 
36000  km  and  inclination  between  20°  and  30°,  with  the  aim  of  avoiding  the  need  for 
numerical  I  :it  egra  t  i  on  with  each  new  launch. 

There  are  four  sources  of  error: 

t  1 )  Lsti.  'ting  the  mass/area  ratio  of  the  satellite; 

Id)  Predicting  the  average  solar  activity  during  the  lifetime; 

t  i)  1'rrors  in  perigee  height  of  the  initial  orbit  used; 

(d)  errors  in  the  estimated  mean  perigee  height,  caused  by  the  approximations. 

Sources  (I)  and  (.2)  would  both  be  expected  to  give  errors  of  10-20%  in  lifetime.  Sources 
(3)  and  (-U  would  both  be  expected  to  give  errors  of  about  10  km  (sd)  in  perigee  height, 
equivalent  to  errors  by  a  factor  of  two  in  lifetime  for  low  perigee  ( 1 30  km  height)  or 
20%  error  for  high  perigee  (250  km).  The  errors  (1)  to  (3)  will  remain,  so  there  seems 
little  point  in  trying  to  avoid  error  (d)  by  lengthy  numerical  integrations.  The  approx¬ 
imate  method  is  to  be  preferred  for  lifetime  predictions  of  newly-launched  objects  in 
low-inclination  transfer  orbits. 
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Lifetime  Lq  of  satellite  of  mass/srea  m/S  =  100  kg/m^  in 
transfer  orbit  of  eccentricity  0.73.  If  m/S  100  ,  multiply 
Lg  by  m/IOOS.  Lunisolar  perturbations  ignored 


Perigee  height  -  km 


Lifetime  Lq  of  satellite  of  maw/area  100  kg/mz  in  transfer 

orbit  with  perigee  height  between  300  and  700  km,  if  the 
average  future  exospheric  temperature  is  800K  or  900K 
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equation  (22) 
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